Non universal diffusion limited aggregation and exact fractal dimensions 
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In analogy to recent results on non-universal roughening in surface growth [Lam and Sander, Phys. 
Rev. Lett. 69, 3338 (1992)], we propose a variant of diffusion-limited aggregation (DLA) in which 
■ the radii of the particles are chosen from a power law distribution. For very broad distributions, 

^\ ' the huge particles dominate and the fractal dimension is calculated exactly using a scaling theory. 

0^ , For narrower distributions, it crosses back to DLA. We simulated 1200 clusters containing up to 

200, 000 particles. The fractal dimensions obtained are in reasonable agreement with our theory. 

<D 
Ph 

Diffusion limited aggregation (DLA) ]lj is an impor- 
tant prototype of fractal growth 0. In two dimensions, 
its fractal dimension is determined to be 1.712 ± 0.003 
from large scale off- lattice simulations . There exists 
• no exact formula for the fractal dimension. In spite of in- 
tensive studies, the task of obtaining accurate estimates 
. of the fractal dimension by analytical means remains to 
be very challenging MM . 

In this work, we investigate a new class of variants with 
exactly solvable fractal dimensions for a certain range of 
a model parameter. In another range, the model crosses 
back to DLA. The transition is associated with an in- 
| teresting morphology change. The model might have 
I . relevance to cluster-cluster aggregation (CCA) models 
'"O ' which simulate colloidal aggregation More impor- 
Ch , tantly, DLA and most previous non-trivial variants are 
not expected to have exactly solvable fractal dimensions. 
\ Our variant might have a particular role in testing the- 
^ ■ oretical approaches of branched growth [^],[|. However, 
kji \ even though the variant has an exact dimension for a cer- 
. tain regime which extrapolates continuously to DLA, our 
C$ ' approach is not capable of calculating the fractal dimen- 
sion, Ddla-, of DLA. The reason is that the transition 
point is not known a priori and has to be expressed in 
terms of Dp la- We concentrate on the two dimensional 
case, while the results can be generalized easily to higher 
dimensions. 

Our variant of DLA is motivated by recent results on 
Zhang's model of surface growth with power-law noise 
Hfr|. The algorithm is very similar to standard DLA. 
Particles are launched one by one and carry out Brow- 
nian motion. Starting from an immobile seed, the ag- 
gregate grows when a walker hits it and becomes part 
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of it. Usually, all the particles launched have the same 
radii. The asymptotic scaling properties and the mor- 
phology are not expected to change if the particle radii 
are random but bounded. However, the universality can 
be broken if the probability distribution of the radii is 
very broad. This leads us to characterize our variant 
by the following power-law distribution, P(r), of particle 
radii, r: 

PM = ( » /r " +1 fOT T > 1 (1) 

[ > \ for r < 1 w 

It recovers the standard DLA algorithm as /i — > oo. In 
general, the walkers can be fractals themselves with di- 
mension Dp. An interesting example is the CCA model 
in which, for some regimes, the radii of the ensemble 
of clusters can have a power-law distribution Qj. The 
walkers in our variant then correspond to the wander- 
ing clusters in the CCA model. In this case, the fractal 
dimension of the walkers and the final aggregate should 
be the same. This property follows naturally from our 
scaling theory in the appropriate regime. 

We focus on the case that the aggregate and the indi- 
vidual particles are dense enough so that they are not 
transparent to each other. This is generally true in 
low dimensions. For programming convenience, we fur- 
ther assume in our simulations that the wandering parti- 
cles have a circular outer boundary. This simplification 
should not alter the scaling properties. 

The simulation of off-lattice DLA with the power law 
distribution of the particle radii in Eq. (|l]) is compli- 
cated due to the absence of an upper cutoff for the parti- 
cle radii. To determine accurately whether a walker has 
touched the cluster, it is no longer sufficient to search in 
a finite neighborhood for centers of other particles. In 
order to obtain an efficient algorithm, the focus must be 
shifted from the centers of the particles to their perime- 
ters, which are instead stored in a hierarchical map struc- 
ture HQ . We can thus check efficiently whether a walker 
with a given radius overlaps with a perimeter site of some 
particle on the cluster. However, this procedure signif- 
icantly increases both the storage requirements and the 
run time of the simulation, especially for a small power 
law exponent \i — ► 1. Therefore, we had to restrict our 
simulations to clusters with masses between 100, 000 and 
200,000 particles. 

Figure |l| shows three typical clusters grown for expo- 
nents fi = 2.5, /j, = 1.713 and /i = 1.3. Each of them 
contains 100, 000 particles. The choice of the values of 
/x is to illustrate a morphology transition which will be 
explained later using a scaling theory. For /i = 2.5 [Fig. 
|l|(a)] , the presence of the large particles results in a clus- 
ter significantly more noisy than DLA. However, the 
small particles are abundant enough to dominate the ge- 
ometry and the overall branching structure is very simi- 
lar to DLA. There exists a number of well defined main 
branches emanating radially from the center of the clus- 
ter. We observe that the ratio between the size of the 
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largest particles and that of the cluster decreases as the 
number of particles, TV, increases. In fact, we will ex- 
plain that the cluster becomes indistinguishable to DLA 
for very large TV. In contrast, for the very broad distribu- 
tion at fj, = 1.3 [Fig. 0(c)], the morphology is completely 
altered by the very big particles. A well defined geomet- 
rical center of the cluster does no longer exist. The size of 
the largest particles is comparable to that of the cluster 
and dominates the structure completely. The morphol- 
ogy is more similar to that of CCA than DLA. We will 
show that the transition between the two different mor- 
phologies occurs at \x = Ddla — 1-713 [Fig. ||(b)]. Here, 
the DLA type morphology holds marginally and the size 
of the largest particles compared to the cluster size de- 
creases very slowly as TV increases. 

We now present a scaling theory for the variant, which 
predicts a transition from a regime dominated by small 
particles to one with dominating large ones as /u de- 
creases. Consider the scaling form 

R G ~ TV 1 / 7 (2) 

where Rq is the radius of gyration of the cluster and TV 
is the number of particles. For cases such as D p — or 
fi —> oo, the exponent 7 reduces to the fractal dimension. 
The general relationship will be worked out later. 

We first examine a broad distribution with \i < Ddla- 
We are going to show that 7 is exactly given by 7 = /i < 
Ddla- The idea of the proof to be presented below can 
be sketched out briefly as follows: If it were true that 
7 > fi, the cluster radius Rq would be even smaller than 
the typical size of the largest particle inside the cluster. 
If 7 < fj,, Rg would be too large to be accounted for 
by the large particles. This leads to a domination of 
small particles, which again can be proved to be wrong. 
Our derivation is closely related to the analogous one for 
surface growth with power-law noise, in which case Lam 
and Sander Q| proved the exactness of a formula for the 
scaling exponent suggested independently by Zhang ||] 
and Krug @. 

We now give our arguments in more detail. It is easy 
to show that 7 ^ Ddla, which necessarily implies a non- 
universal behavior fj, < Ddla- For a cluster of TV par- 
ticles, the expected radius of the largest particle, r rnax , 
sampled from the power-law distribution in Eq. (|l|) , fol- 
lows r max ~ TV 1 '' 1 . The proposition 7 = Ddla implies 
Rq ~ n 1 / Ddla . This is impossible for fi < D D la be- 
cause at sufficiently large TV, the largest particle would 
even be bigger than the cluster itself (r max > Rg)\ The 
same argument shows more generally that 7 > /i is false 
when < Ddla- 

In fact, the huge particles are not only relevant to the 
scaling, as proved above, but also dominate the geome- 
try completely. It is plausible that the largest particle 
sets the scale of the cluster: Rq ~ r max . This assump- 
tion immediately leads to 7 = fj,, since r rnax ~ N 1 ^. 
We will show that this simple picture is indeed correct 
asymptotically. 
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The conventional way to compare the visual appear- 
ance of DLA of different sizes is to rescale them to a 
standardized radius of gyration. Instead, we now rescale 
the clusters by a factor N^ 1 ^ and examine the depen- 
dence of the radii of the clusters on JV. In this rescaling 
scheme, the size of the biggest particle is independent of 
N. Furthermore, the radii, r', of the rescaled constituent 
particles follow a new probability distribution P'(r'). In 
a cluster of N particles, the number density of particles 
per unit r' is n(r') = NP'(r') given by 

where the lower cutoff radius, r m , is a function of N given 
by r m (N) = N -1 '^. This implies that not only r max but 
also the whole distribution n(r'), except for the lower 
cutoff r rn , is independent of N. The simple picture of 
the domination of the large particles corresponds to the 
assumption that the variations of the cutoff can be ne- 
glected asymptotically. If this is true, n(r') is completely 
independent of N as N — > oo and both, the self-similarity 
and 7 = fx, become obvious. 

For narrower distributions with /i > Ddla, the small 
particles dominate. We have 7 = Ddla and the mor- 
phology is the same as DLA. The huge particles are irrel- 
evant to the geometry. The radius of the largest particle 
vanishes compared with the cluster size as N increases. 
In summary, we have: 

7 = min{/x, D DLA } (4) 

Although the complication about the lower cutoff does 
not alter the exponent, it leads to strong crossover ef- 
fects. Using the above rescaling scheme with the factor 
-/V -1 /^, as N — > 00 the rescaled cluster radius, R' , con- 
verges for ,u < Ddla but diverges for \i > Ddla- The 
H = Ddla case is marginal and R' is expected to diverge 
logarithmically, and by analogy with the results for sur- 
face growth we postulate a logarithmic correction to 
scaling at /i = Ddla- 

Rg ~ N 1/>Dr ' LA [log(N)} 1 / 2 (5) 

Now, we compute the fractal dimension D of the clus- 
ter from the relation M ~ Rq, where M is the mass 
of the aggregate. Using the distribution of the particle 
radii in Eq. (|lj) and the assumption that the constituent 
particles are fractals of dimension D p , we get 

M ~ J N D p < fi 

\ Lt, < Dp (6) 

where the averaged particle mass diverges for small fi. 
At fi = Dp, the divergence is marginal: 

M~NlnN (7) 

Combining Eqs. (|2|), (|4|) and (|^), the fractal dimensions 
are obtained for four different cases: 
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Ddla D p ,D D la < fi 

D = ) D p D DLA /fi D DLA < fi < D p , , 

P D p < n < D DLA { > 

Dp < D p ,D DLA 

The variant crosses back to DLA for narrower distribu- 
tions of Dp, Ddla < M- F° r Ddla < M < -Dp, even 
though the morphology is the same as DLA, the fractal 
dimension is different since the averaged particle mass 
diverges. For the other two cases of broad distributions, 
the morphology is non-universal. It is particularly in- 
teresting that for \i < Dp, Ddla, D = D p follows. It is 
precisely how it should be when the walker represents the 
clusters in the CCA model. This sets an upper bound 
for /i if the variant describes some regime of CCA. Sim- 
ulations of CCA in general give n or its effective value 
well within this bound Q . 

We test our scaling theory numerically for the case 
Dp — 2. Equation (||) then reduces to 

2 < fj, 

Ddla < M < 2 (9) 
/i < Ddla 

The verification of Eq. (^) also establishes the validity 
of Eq. (f|) and the general result in Eq. (|J) , since they 
are all related by the trivial Eq. (||). The logarithmic 
corrections in Eqs. (||) and (0) at the transition points 
give respectively: 




M~l Rg DLA /Io9(N) v^Ddla M()1 



In our naive computation of D from the slope of the best 
fitted straight line in the log-log plot of M against Rq, 
these corrections cause an underestimation at fi = Dd la 
and an overestimation at /i = 2. 

Let Ti and Xi be respectively the radius and the po- 
sition of the i-th particle in a cluster with N particles. 
Assuming that every particle has uniform unit density, 
the mass of the individual particle is m, = irrf. Since 
the particles are disks instead of points, the mass and 
radius of gyration of the cluster is given by: 

N 

M = m i (H) 

i=l 

N A 
MRl=J2 m ^~^M) 2 + ^} (12) 



where xqm is the center of mass of the cluster. 

We have grown a total of 1200 clusters for exponents 
/i in the range 1.0 < /x < 5.0. For various cluster sizes 
N, we compute the ensemble averaged cluster mass and 
radius of gyration. Each data point in the Rq vs M plot 
is obtained from ensemble average over clusters of fixed 
N. However, the averaging has to be done with caution 
to avoid divergence. For example, Eq. (Q) implies that 
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the probability distribution, P(m), of the mass, m, of 
the walkers follows the power-law: P(m) ~ 1/?tW 2+:l . 
For [i < 2, the arithmetic mean of m diverges. However, 
the geometrical ensemble average, exp < InM > is well 
defined, where the bracket denotes the arithmetic mean. 
Similarly, we also take the geometrical average exp < 
In Rg > for the radius of gyration Rg ■ Figure (||) shows 
the geometrical ensemble averages M against Rq in a 
log-log plot for the selected values \x = 1.3, \x = 1.7, and 
/i = 2.5. The numbers of clusters used are 50, 104 and 
25 respectively. For all values of /z we investigated, we 
obtain a reasonable scaling behavior. 

Figure ^| shows the measured fractal dimensions D as 
a function of /i. We computed D by averaging over the 
dimensions for each individual cluster, obtained from the 
corresponding scaling plot of Rg against M. The er- 
ror bars were obtained from the statistical fluctuations. 
Quantitatively the same result is obtained when we com- 
pute D from the slopes of the geometrically averaged Rg 
vs M plots (Fig. ||). Also shown in Fig. || is the predic- 
tion of the scaling argument in Eq. @ . Good agreement 
with our theory is observed far away from the transition 
points for /i ~ 1 and fi > 2.5. At the transition points 
/i = Ddla and 2, the expected discrepancies due to the 
logarithmic corrections in Eq. ( |l0| ) are observed. 

In summary, we propose a variant of DLA in which the 
random walkers have random radii of very broad power 
law distribution. A scaling theory predicts that, as the 
distribution becomes narrower, there exists a transition 
between regimes dominated by large and small particles. 
The fractal dimension can be calculated exactly for the 
former regime. At the transition point we find a logarith- 
mic correction to scaling. An analogous scaling theory 
leading to a similar transition has been verified numeri- 
cally for the case of surface growth with power law noise. 
We generated 1200 clusters of our DLA variant. When 
huge particles dominate, the morphology of the cluster 
is different from that of DLA and might be related to 
that of the cluster-cluster aggregation model. The mea- 
sured fractal dimension is in reasonable agreement with 
the scaling theory. The numerical precision is limited 
by strong statistical fluctuations and finite size effects, 
which manifest themselves as logarithmic corrections at 
the transition points. 
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FIG. 1. Typical clusters at (a) the small particles domi- 
nating regime (/i = 2.5), (b) the transition point (/i = 1.713), 
and (c) the large particles dominating regime (fi — 1.3). 

FIG. 2. Dependence of the cluster mass M on the radius 
of gyration Rq for three selected values of ji. The geometrical 
mean is adopted for the ensemble averaging. 

FIG. 3. Dependence of the fractal dimension D on the 
exponent fx. The solid line denotes our numerical results, 
whereas the dashed line marks the prediction of the scaling 
theory. The systematic discrepancies at the transition points 
fi ~ 1.713 and fi = 2 are due to finite size effects in the form 
of logarithmic corrections. 
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